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A PROCEDURE FOR RAPID COMPUTATION OF EQUILIBRIUM 
CHEMICAL COMPOSITION AND THERMODYNAMIC PROPERTIES OF THE 
GENERALIZED ARGON-NITROGEN-OXYGEN-HYDROGEN GAS SYSTEM 

By Vincent R. Mascitti 
Langley Research Center 

SUMMARY 

A procedure is described for the rapid computation of the equilibrium chemical 
composition and thermodynamic properties of the generalized argon-nitrogen-oxygen- 
hydrogen gas mixture system for nonionized conditions. The procedure employs a varia- 
tion of the Newton iteration scheme to solve for chemical composition. The generalized 
model reduces to the hydrogen-air, hydrogen-oxygen, pure hydrogen, and equilibrium air 
systems. A computer program listing of the calculation procedure for this gas model is 
presented. 


INTRODUCTION 

The calculation of thermodynamic properties of combustion gas mixtures at high 
temperatures is complex, difficult, and time consuming. These difficulties arise pri- 
marily because of the determination of equilibrium chemical composition which requires 
the solution of a system of equations that are not simultaneously linear. In the past, iter- 
ative methods of solution such as those summarized in reference 1 have been used exten- 
sively. However, these methods have proved to be so time consuming as to preclude their 
direct use in calculations of engine performance, flow -field properties, and kinetic pro- 
cesses. As a result, the thermodynamic properties of combustion gases calculated by 
these methods have been presented in tabular form in the literature. For example, refer- 
ence 2 is an elaborate treatment of the hydrogen-air system. To apply the results of ref- 
erence 2 to more complicated problems requiring solutions by computer, the user is com- 
pelled to fit curves to the gas mixture thermodynamic properties; thus, this technique 
sacrifices generality, computer storage, and, most important, accuracy. 

In reference 3, a new method was developed to calculate chemical composition. 

The method employs the Newton iteration scheme and is as much as two orders of mag- 
nitude faster in computational time than previous methods. A variation of the Newton 
iteration scheme has been applied with success to the simplified hydrocarbon-air com- 
bustion gas model of reference 4. The assumption of reference 4 is that atomic nitrogen 



species form in negligible amounts in the range of operation of air-breathing propulsion 
engines and, therefore, have a small effect on the performance of these engines. How- 
ever, the formation of atomic nitrogen species, such as ON and N, is important for hyper- 
sonic engine nozzle calculations which include the effects of chemical kinetics (ref. 5). 
During the study of reference 4 it was found that the atomic nitrogen species could be 
included within the framework of the same iteration scheme if carbon species were 
neglected. 

The considerable improvement in speed of solution realized by the use of the Newton 
iteration scheme for the simplified hydrocarbon-air combustion gas model in reference 4 
makes it desirable to extend the procedure to other gas models. The gas model of the 
present study includes the formation of the following chemical species: 

O, N, Ar, H, 0 2 , N 2 , H 2 , ON, OH, H 2 0 

The computer program presented in appendix A is geared toward a generalized approach 
to the argon-nitrogen-oxygen-hydrogen system in that the initial proportions of atoms 
are input variables. This procedure allows the automatic reduction of the general gas 
model to the hydrogen-air, hydrogen-oxygen, pure oxygen, pure hydrogen, and equilibrium 
air systems. This feature permits the effect of secondary injection of oxygen or hydrogen 
peroxide to be taken into account. Thrust augmentation by secondary injection has been 
considered for some hypersonic engine cycles. 

The computer program can be readily incorporated as a subroutine to a larger pro- 
gram or used alone to generate Mollier diagrams. Inputs to the program are the ratio of 
nitrogen atoms to oxygen atoms, the ratio of argon atoms to oxygen atoms, equivalence 
ratio, temperature, and pressure. 


SYMBOLS 

A,B coefficients of iteration functions 

C = 1 + 2r^ r + r N 

Cp specific heat at constant pressure, Btu/lbm-°R (J/kg-°K) 

G Gibbs free energy, Btu/lbm (JAg) 

H enthalpy, Btu/lbm (J/kg) 
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t 

K equilibrium constant in terms of partial pressures; numerical subscripts 

indicate individual reaction 

M mean molecular weight 

p pressure, atm 

R universal gas constant, 1.98588 Btu/mole-°R (8.314 J/mole-°K) 

r Ar ratio of argon atoms to oxygen atoms 

r e q equivalence ratio 

r ^ ratio of nitrogen atoms to oxygen atoms 

S entropy, Btu/lbm-°R (J/kg-°K) 

T temperature, °R (°K) 

X mole fraction 


Subscripts: 

I Ith temperature 

J Jth species 

Superscript: 

o at reference pressure of 1 atmosphere (1.01325 X 10^ N/m^) 

ANALYTICAL APPROACH 

The general equations presented in this analysis for computing chemical composition 
and combustion gas properties make use of the following assumptions: 

(1) The derived gas mixtures simulate the products of complete, constant-pressure, 
adiabatic combustion. 

(2) The gas mixtures are in thermal and chemical equilibrium. 
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(3) The ideal gas equation of state is valid and only gas phases are considered. 

(4) The effects of ionization are ignored. 

With these general assumptions the determination of chemical composition must be 
accomplished in the presence or absence of dissociation. Once the chemical composition 
is determined and the thermodynamic properties of the pure constituent are provided, the 
thermodynamic properties of the gas mixture can be established. 

Chemical Composition 

The chemical species assumed for the gas model of this analysis are 

O, N, Ar, H, 0 2 , N 2 , H 2 , ON, OH, H 2 0 

The solution of this system is centered about the Newton iteration scheme as suggested 
in reference 3. However, the inordinate amount of labor required to reduce the equations 
which represent this system to a single iteration equation was impractical. As a result, 
the equations were reduced to two equations and two unknowns. Newton's iteration scheme 
is then applied to the two equations simultaneously, in a manner similar to that presented 
in reference 4. It should be emphasized that rapid convergence is not necessarily assured 
when employing the Newton iteration scheme becuase of the highly nonlinear character of 
the equations. With this observation in mind, the governing equations for the solution of 
chemical composition can now be developed. 

In a system of 10 chemical species and four chemical elements, there are six 
independent chemical reactions which can be written with the corresponding equilibrium 
expressions as 


20 = 0 2 

(1) 

Xq 2 = P K 5 Xq 

2N= N 2 

(2) 

X N2 = pKgxjsf 

2H = H 2 

(3) 

X H2 = PK 7 X| 

O + N = ON 

(4) 

X ON = P k 8 X 0 X N 

O + H = OH 

(5) 

x OH = P K 9 X 0 X H 

O + 2H = H 2 0 

(6) 

Xh 2 0 = P 2 KioXqXh 
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Application of the conservation of mass for each chemical element gives 

(7) £x J= l 
J 

(8) X N + 2X N2 + Xqn = r N( x O + 2x 02 + X 0N + x OH + x H 2 o) 

( 9 ) X^ = r^JXQ + 2 Xq 2 + X ON + X OH + X H2 o) 

(10) X H + 2Xjj 2 + Xq H + 2Xjj 2 q = 2r e q^X 0 + 2X 0 2 + X ON + x OH + X H 2 o) 

Solving for two equations and two unknowns in Xq and Xjj yields 
a = A 4 0 Xq + (a 32 x| + A 31 X h + A 30 )x^ + (a 24 x£ + A 23 X 4 + A 22 x| + A 2 iX h + A 2 o)x% 

+ (ai 4 Xh + A 13 x|[ + A 12 x| + AhX h +Aio)xq + (aq 4 Xh + Aq3Xh +A()2Xh + AoiXh + Aoo) = 0 

/3 = ^B 22 Xq + B 21 Xq + B 2 ojXjj + ^Bi 2 Xo + Bi ]Xq + Bio)xh + (bo2Xo + BoiXo) = 0 
where 

A 40 = C 2 p3K 5 (4K 5 K 6 - Kg) 

A 32 = Cp 4 K 10 |4(l + C)K 5 K 6 - (l + rAr )K| 

A 31 = Cp3K 9 [4(l + C)K 5 K6 - (1 + r^xf] 

A 30 = Cp 2 [ 4(1 + C)K 5 K6 - (1 + r^Kj + 2( rAr - C)K 5 K 8 
A 24 = (1 + C) 2 p5K6K? 0 
A 23 = 2(1 + C) 2 p4K 6 K 9 K 10 

A 22 = p3 ^(1 + C) 2 K 6 K 1q + 8CK 5 K 6 K 7 + (1 + C) 2 K 6 K 9 + C |c - 2(l + + K 7 k1 

+ [2C(l + rAr ) - (1 + C)(l + C + r^KgKioj 


5 


A 2 i = P 2 jsCKgKg + 2(1 + C) 2 KgKg + c|<t - 2(l + 

+ J2C^1 + r ^ - (1 + C)(l + C + r Ar)J K 8 K 9| 

a 20 = p 2 (c|2(i + rAr ) - c]k| - 8CK 5 K 6 j 

+ p^l + 2r M - 2c)k 5 + (1 + C) 2 K 6 + 2C(1 + r^) - (1 + C)(l + C + rAr) K 

A 14 = 4(1 + C)p 4 K 6 K 7 K 10 
A13 = 4(1 + C)p 3 K6 (k 10 + K 7 K 9 ) 


Ai 2 = p 2 


, 4{1 + CjKgK-j + 4(1 + C)K 6 K 9 + ( rAr - c)k 10 - 2(1 + rAr )K 7 K 8 ] 
-4p3(l + C)K 6 K 1q 


An = P 4(1 + C)Kg - 2^1 + r^Kg + - c)Kgj - 4(1 + C)p 2 KgKg 

A10 = P J2 (l + r^Kg - 4(1 + C)Kgj + (r^ - c) 

A 0 4 = 4p 3 KgK 7 

A03 = 8p 2 K 6 K 7 

A 02 = p(4Kg - K ? ) - 8p 2 K6K 7 

Aqi = -1 - 8pK 6 


Aqq = 4 P k 6 + 1 
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and 


B 2 2 = " 2 ( r eq + c )p 3k 8 K 10 

B 21 = -2p 2 |(2r eq + C)K 7 K 8 + (1 - r eq )K 10 ] 

B 2 o = -2pK 7 

B 12 = “( 2r eq + C)P 2 K 8 K 9 

B 11 = -(l - 2r eq jpKg - ^4r eq + C^pK 8 

B 10 = -i 

b 02 = - 2r eqP( K 8 ~ 2K 5) 

B 01 = 2r eq( 2 P K 8 + x ) 

C = 1 + 2r Ap + r ^ 

Since the expression for /3 is second order, it can be satisfied by the quadratic 
formula if one unknown is assumed. The resulting pair of values ^Xq,Xjj^ can then be 
tested for compliance with a within the desired accuracy. Although there may be a 
number of pairs of values ^Xq,Xjj) which satisfy a = 0, only one pair satisfies the condi- 
tion 0 S Xj s 1. For r eq S 1, Xq was chosen to be the iteration variable. 

Newton's iteration scheme in two variables can be written as 

Xo(L t „ . x^L, - f ° (L) ’ XH - ( ^ , 

j!-(*[xo<L),Xh<L)]J 

where Xq(L) is the Lth approximation of Xq and the derivative of a^Xo,Xjj) with 
respect to Xq is 
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d 

dX Q 


a 


X g (L),X h (L) 


da 

dXr, 



/»/3 A 

da 

f 3Xo\ 

3X h 

1 9/3 


\3Xti / 


For r g q > 1, Xjj was chosen to be the iteration variable since a and /3 were insen- 
sitive to changes in Xq in this region. 

Because of the highly nonlinear character of this system and the presence of extra- 
neous solutions in the neighborhood of the correct solution, an arbitrary choice of an ini- 
tial value for the iteration scheme was wholly inadequate. Rapid convergence by Newton’s 
iteration scheme can be realized only when the function is nearly linear or an initial esti- 
mate is chosen in the neighborhood of the solution. As a result, considerable intelligence 
has been built into the selection of the initial value of the iteration to insure that it is 
approximately equal to the correct solution over the entire range of variables. 

The choice of an initial value in the neighborhood of the solution has been made by 
formulating a simplified gas model for a particular region of equivalence ratio. The sim- 
plified models chosen were those which would provide a solution in a closed-form expres- 
sion. For example, in the region r e q < 1, the following simplified gas model was chosen: 


h 2°> °2> °> n 2> ^ 


In general, the simplified models chosen for the initial iteration values give results 
nearly equal to the solution for the entire system. However, correction factors were nec- 
essary to insure convergence to the correct solution. The simplified gas models and cor- 
rection factors are developed in detail in appendix B. 

Below T cu t_ 0 ff (T cu t_off = 100 logjo P + 2800), the gas mixture is assumed to be 
nondissociating. The cut-off temperature is an arbitrarily defined temperature below 
which the effects of dissociation can be neglected. Below this temperature limit, the ini- 
tial proportions of elements are sufficient to define the gas composition. 

Chemical species assumed for the simplified gas model in the region T 2i T cut _ Q fj 

are 

A r> °2> n 2> h 2> h 2° 

and the mass conservation equations are 
J 

2Xn 2 = r N (2X 02 + Xh 2 o) 
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''Ar “ r Ar 


pXo 2 + Xh 2 o) 


2Xh 2 + 2x h 2 0 “ 2r eq( 2x 02 + Xh 2^) 


Although there seems to be one less equation than required, in each region of equivalence 
ratio one or more species can be neglected. For r gq i 1, 



x h 2 o- 


2 r eq 
c + r eq 


For r eq > 1, 


X O2 = 0 


X 


_ 2 
h 2° C + 2r eq 


- 1 


At low temperatures for fuel -rich mixtures ammonia may form in substantial 
amounts. In the present study, the formation of ammonia is handled just as is done in 
reference 4, that is, the computer program indicates the formation of ammonia by an 
error statement. When the limiting pressure for ammonia formation is exceeded, the 
error statement is printed, but the program computes chemical composition as though 
ammonia had not formed. 


Thermodynamic Properties of Pure Constituents 

The standard reference states of the elements Ar, N, O, H taken from reference 2 
are as follows: 

Ar as Ar 
N as N 2 
O as O 2 
H as H 2 

By definition, the energy content of these elements in their standard reference states 
(T = 0° R) is zero. 


The specific heat at constant pressure for each species 


: P,J 

R 


is also taken from 


reference 2 and tabulated in the computer program as a function of temperature. The 

thermodynamic properties of the pure constituents U— -^-1 are obtained by numerical 
integration in the form 


R 


R 


'P- 


I.J 


+ c r 




I-l) 
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J - tf-1, J , 

R R 



The thermodynamic properties of the pure constituents obtained in this manner agree with 
those tabulated in reference 2 within four significant figures. With the thermodynamic 
properties of the pure constituents defined, the equilibrium constants can be formed. For 
a general reaction, the equilibrium expression is 


aA + bB = cC + dD 


and the equilibrium constant is 


K = exp 


c + d £e 

RT RT 




where 


Gj Hj $ 
RT RT " R 


Thermodynamic Properties of Gas Mixtures 

The foregoing equations have outlined the solution for gas mixture composition. 

With the chemical composition defined and the thermodynamic properties of the pure con- 
stituents provided, the thermodynamic properties of a gas mixture can be calculated by 
the following equations (ref. 2): 


M = 


I 


X jMj 




where 



5 

R 


are the thermodynamic properties of the pure constituents at refer- 


ence pressure (1 atm). These thermodynamic expressions inherently assume that the 
pure constituents obey the perfect gas law. 
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A listing of the computer program to calculate the composition and thermodynamic 
properties of the argon-nitrogen-oxygen-hydrogen gas model is presented and discussed 
in appendix A. Computational time for obtaining the chemical composition and thermody- 
namic properties with this program has been estimated at 12 000 cases/min on the Control 
Data 6600 computer system at the Langley Research Center. 

RESULTS 

The computer program for the present gas model was run over a wide range of input 
variables to establish its reliability. All combinations of the following variables were 
computed: 


Variable 

Range 

Increment 

r N 

0 to 5 

1 

r Ar 

0 

0 

r eq 

0 to 1 

0.1 

1 to 9 

1 

T 

T cut-off to 10 000° R 

200° R 


(5556° K) 

(111° K) 

P 

0.001 to 100 atm 

Orders of magnitude 


For all these variable combinations, the iteration procedure converged to the proper solu- 
tion with an over -all average of five iterations per case. 

This program has been used to calculate the chemical composition and thermody- 
namic properties of hydrogen-air combustion. The results are presented in tables I, II, 
and HI for equivalence ratios of 0 (air), 1, and 2. 5, respectively. Ranges of temperature 
and pressure are also included. Results from the tabulations of reference 2 are presented 
for comparison. The good agreement between gas models is not surprising, since the 
present model assumes the same chemical species with the exception of NH and NHg, both 
of which form in negligible amounts. 

The values in table m for r eq = 2.5, T = 6000° R (3333° K), p = 0.001 atm and 
r eq = 2.5, T = 10 000° R (5556° K), p = 100 atm show substantial disagreement with 
those of reference 2. In a localized temperature range for fuel -rich mixtures at extreme 
pressures, the value of Xq in the quadratic expression for /3 is imaginary. For these 
cases, the radical term is deleted so that solutions can be obtained. The computer pro- 
gram prints an error statement to alert the user that results should be viewed with caution 
since they may be substantially in error. 



CONCLUDING REMARKS 


A rapid computational procedure has been developed to determine the equilibrium 
chemical composition and thermodynamic properties of the generalized argon-nitrogen- 
oxygen-hydrogen gas system. This system includes the formation of the following chem- 
ical species: 


O, N, Ar, H, 0 2 , N 2 , H 2 , ON, OH, HgO 

Newton's iteration scheme is employed to calculate chemical composition with initial 
estimates being obtained from the solution of a simplified gas model. Chemical compo- 
sition is determined within five significant figures with an average of five iterations being 
required per case. 

The reliability of the computational procedure was established over a wide range 
of input variables. The computer program for this gas model has been used to calculate 
the chemical composition and thermodynamic properties of the hydrogen-air system for 
a range of variables. The results agree with those obtained by more comprehensive 
treatments. 

The automatic reduction of the generalized gas model to such systems as hydrogen- 
air, hydrogen-oxygen, oxygen, hydrogen, and air facilitates the computation of thermody- 
namic processes involving the secondary ejection of additional fuel and/or oxidizers. 

This feature should prove useful when determining the performance of air-breathing engine 
cycles which employ forms of thrust augmentation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 12, 1969, 

720-03-10-02-23. 
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APPENDIX A 


COMPUTER PROGRAM FOR CALCULATING CHEMICAL COMPOSITION AND 
THERMODYNAMIC PROPERTIES OF THE Ar-N-O-H GAS SYSTEM 

The calculation procedure described in the main body of the paper for obtaining the 
chemical composition and thermodynamic properties of this reacting gas mixture has been 
programed for high-speed digital computation. The purpose of this appendix is to provide 
a description of the necessary input and available output as well as a FORTRAN IV (ref. 6) 
listing of the source program. An example input case and the resulting output listing are 
included. 


Description of Program 

The program reads in the necessary reference enthalpy and Gibbs free energy and 
the specific heat at constant pressure of each chemical specie as a function of tempera- 
ture. The variation of enthalpy and entropy for each specie at reference pressure (1 atm) 
is generated over the full temperature range by numerical integration. The program then 
reads the initial proportions of atoms, temperature, and pressure. With temperature 
known, the thermodynamic properties of the individual species are interpolated linearly. 
The Gibbs free energy, equilibrium constants, and the coefficients of the iteration func- 
tions are computed. The iteration for chemical composition is performed. With the 
chemical composition defined, the thermodynamic properties of the gas mixture are then 
calculated. 


Program Listing 

The FORTRAN listing of the source program used at the Langley Research Center 
on the Control Data series 6000 computer systems is as follows: 
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C THERMODYNAMIC PROPERTIES OF ARGON-MI TROGEN-OXYGh N-HYDROGEN SYSTh ■n 
CO N AP H O 2 N2 H2 OM OH H20 

C X < 1 ) X ( 2 ) X. ( 3 ) X ( 4 ) X(5> X (6) X ( 7 ) X(8) X<9) X(10) 

DIMENSION X (Q9) ,Y (99) .7 (99) ,F (99) .CP ( 1 10. 10 ) »H ( 1 10 . 10 ) .HT ( 1 1 0 . 1 0 ) . 
1 S ( 1 1 0 » 1 0 ) . ST ( 1 1 0 . 1 0 ) . G ( 1 1 0 , 1 0 ) » GT ( 1 1 0 . 1 0 ) . T T (11 0 ) . HR ( 1 0 ) . 

2SR (10) . GRT (10). MR (10) 

DATA (MR (J) »J= 1.101/16.. 14. 00B. 19.944. 1 . 008 , 32. . ?H . 0 lft . 2. 0 1 ft . 

1 30. 00B. 17. 008. 18. 01 ft/ 

REAL K5.K6.K7 .K8.K9.K1 0 . KM , MR 
NAMELIST/NUM1 /HT.GT.CP 
N A MEL I ST/NUM2/RN .RAR.RFO.T .P 
READ (5.NUM1 ) 

J=1 

43 1 = 1 

TT < 1 » =200. 

H ( 1 , J) =HT ( 1 . J) *200. 

G ( 1 » J) =GT ( 1 . ,)) 

S ( 1 . J) =HT ( 1 . J) -GT ( 1 . J) 

ST < 1 . J) =S ( 1 . J) 

41 IF ( I. EO. 10ft) GO TO 42 
1 = 1 + 1 

TT ( I ) =TT ( T-l ) +1 00. 

H ( I • J ) =H ( I - 1 . .) ) + ( CP ( I . J) +CP ( l - 1 . . J ) ) /2 . * ( T ( (I)-Tl ( I — 1 ) ) 

S( I » J) =S ( T-l . J) + (CP( I . J) +CP( r-1 , J) ) * ( TT ( I ) - 1 r ( I -1 ) ) / ( 1 T ( I > + TT ( ! - I ) 
1) 

HT ( I . J) =H ( I , J) /TT ( T ) 

ST (I* J)=S( I. J) 

GT ( I . J) =HT ( I . J) -ST ( I . J) 

GO TO 41 

42 IF ( J.EU. 1 0 ) GO 10 SOO 
J=J+1 

GO TO 43 
500 READ (5.NMM2) 

100 FORMAT (SFlft. 8) 

IF (EOF. 5)400. 401 

400 STOP 

401 DO 12 1=1 .12 

12 X ( J)=0. 

PRINT 300.RN.PAR.RFQ. 1 .P 

30 0 FORMAT (8X3HRM=F 15. R. 4X4 HP AR = F 1 -> . 8 . 4 X4HRF0=I- 1 5 . m/2 0 < 2h T=h J 5 . s . 1 b R . 
14X2HP=F 15.B.3HATM/) 

DO 44 1=1.110 
IF (T-TT ( I ) ) 48.45. 44 

44 CONTINUE 

45 K= I — 1 

DFLH= (T-TT (K) )/(TKI)-TT(K)) 

DELS= ( A[.OG ( T/TT (K) i ) / ( A LOG ( TT ( 1 ) / T1 ( K ) ) ) 

DO 46 1=1 .10 

HR ( J) = (HT (K. J) + (H r ( I . J) -HI (K. J) ) *nFI.H) * I 
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SR(J)=ST (Kt.l) + (ST (1 « . I ) — S T ( K » J ) ) *0F LG 
GRT < J) =HP ( J) /T-SR < J) 

46 CONTINUE 

C=1.+2.*PAR+RN 

IF (T .LT. ( 100 . *ALOGl 0 (P) +2«O0 . ) ) GO TO 30 
IF (REO.EO.O.)GO TO 1 
K7=EXP (GRT (4) *2. -GPT (7) ) 

K9=EXP(GRT ( 1 ) +GRT (4) -GPT (9) ) 

K10=EXP (GRT (4) *2. + GRT ( 1 ) -OR I ( 1 <■ > ) 

GO TO 2 

1 K7=0. 

K9=0. 

K I 0 = 0 . 

2 IF (RN.Fq.O. )GO TO 3 
K 6=EXP (GPT (2) *2. -GRT (6) ) 

K8=EXP ( GRT < 1 ) +GRT ( 2) -GRT < 8) ) 

GO TO 4 

3 K6=0. 

K8 = 0 . 

I F (RFO.GT. 20 . ) X (4) = ( SORT ( 1 . +4 . *0**7 ) -1 . ) / ( 2 . *t>*i' / ) 

IF (RFO.GT. 20. )X(7)=1.-X(4) 

IF (RFO.GT. 20 . ) GO TO 11 

4 KS = EXP (GPT ( 1 ) *2. -GRT (S) ) 

A40=C*C*P** 3*«S* (4 . ttXh-K 8*KH/K~> ) / I 0 . **20 

A32=C*P**4*K1 0* (4 • * ( 1 . + C ) * K 6 — ( I . + PAR ) *K8*K8/KS ) / I u .**20 
A3] =C* D **3*K9* ( 4 . * ( 1 . + C ) * K6— ( 1 • + IJ AP ) * KH *KH /KS ) / 1 *29 
A 3 0 = C * P * * 2 * ( 4 . » ( 1 . + C ) *K 6— ( 1 . + R AP ) *KH*K8/KS+2.* ( A R-C ) ) /I 

A24= < 1 . + C ) **2* P **S*K6**< 1 0*K 1 0 / ( KS* 1 0 . **20 > 

A 2 3= 2 ] . + C ) **2* P **4*K6*K9*K ] 0 / ( KS* 1 0 . **20 ) 

A22=P**3* ( 2 . * ( ] . + C ) **2*K6*K 1 ft / ( KG* I 0 . ** 20 ) + R . *1 *a 6*K7 / I 0 . ** 21) + ( l . 

\ +c ) **2*K6*K9**9/ ( KS*l 0 . **20 ) + C* ( C-2 . * ( l . + R AP ) ) *«7* K8 *k 8/ ( KS* 1 o . 
2**20) + (2.*C* n . + R A P ) - ( 1 . +C) * ( 1 . +C + RAR) ) *K M*K 1 0 / ( K 6* I 0 . **2o ) ) 

A 2 1 =P**2* ( 8 . *C*K6/ 1 0 . **20 + 2 . * ( l . + C ) **2*Kb*K9/ ( * b* 1 0 • **2 0 ) + C * ( C- 2 . * 

l(l.-*-RAR))*KR*K8/(KS*10.**20) + (2.*C*(].+RAR)-(l.+C)*(l.+(..+RAR)!*<s 

2 * K 9 / ( K S * 1 0 . * * 2 0 ) ) 

A20 = (P**2* ( C* < 2 . * ( I . + RAR ) -C.) *K8*KR/K C >-K . *C*K6 ) + P* ( ( 1 . + 2 . *RAR-2 . *C ) 

1 + ( ] . + C ) * * 2 * K 6 / K S + < 2 . * C * (1 . + R A R > - ( 1 . ♦ C ) * ( 1 . + C ♦ R A R ) ) * K 8 / K S ) ) / I 0 . * * 2 1 1 
A 14 = 4 . * ( 1 . +C) *P**4*K6*K7*K] 0/ ( KS*1 0 .**20 ) 

A 1 3=4. * ( 1 . + C ) *P**3*K6* ( K 1 0 + K 7*K9 ) / ( KS* 1 0 . * *20 ) 

A]2=(P**2*(4.* (l.+C) *K6*K f/Kb+4 . * ( 1 . +C ) *K6*K.9/K b + ( R A R-C ) *K 1 O/KS-2 . 
1* ( 1 . +RAR) *K7*K8/KS ) -4 . *p**3* ( 1 . + C) *K6*K 1 O/KS ) / 1 0 . **20 
A11=(P*(4.*(1 ,+C)*K6/KS-2.*( 1 . + RAR)*K8/KS + ( R A R-C ) *K9/Ks ) -4 . * ( 1 , *(' ) 
l*p**2*K'6*K9/KS ) /I 0 .**20 

A10=(P*(2.*(1 .+RAR)*K8-4.*(1.+C)*K6) + (RAR-C) ) / ( Kb* 1 0 • **2o ) 
A04=4.*P**3*K6*K7*K7/(Kb*l 0.**20) 

A 0 3 = R . * P * * 2 * K 6 * l< 7 / ( KS* 1 0 . **20 ) 

A02= (P* (4.*K6-K7) /K5-8.*P**2*Ko*K7/KS) /10.**2o 
AO 1=-1 .* ( 1 . +R.*P*K 6) / (K5*1 0 . * * 2 0 ) 

AO 0= ( 4 . *P*K6 + 1 .) / ( K S * 1 0 • * * 2 0 ) 

822 = 2.* (REO + C) * p * * 3 * K 8 * K 1 0 / K S 

(32 1=2. *P**2* ( ( 2.*REO + C) *K7*K8/KS+ ( 1 . — REO ) *K 1 0 /k S ) 
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820=2. «P*K7/K5 

B12=(2.*PE0+C) »p*»2*Ka*K9/K5 

B11=P» ( (l.-2.»REQ> *K9 + ( 4. *RFQ+C ) *K8) /KB 

810=1. /K5 

802=2. ^REQ^P* (K8/K5-2. ) 

B01=-2.*REQ*(2.*P*K8+1 . >/K5 
L=1 

IF ( RED- 1 • ) 6 • 5 « 1 .3 

6 V (1 )= (SQRT ( ( 1 . +C) **2+ l6.*P*Kb* ( 1 . -REO) * (C+REO) ) - (1 ,+C) ) 

1/ (4.*P*K5* (C+REO) ) * ( 1 . 1-.5E-4* I ) 

GO TO 7 

5 Y<1)=1. /(<?.*(! ,+C) *P**2»K1 0) ** U ./ 1.) ) * ( .34-.?K-4*T> 

GO TO 7 

13 Zll) = (SORT ( ( 4 . *REQ + 1 . + C ) «#? +32 . *PFO» ( 2 . *REO + 1 . + C ) *P*K7 ) - ( 4 . *RFO + 1 . 
1+0 ) / <4.*P*K7* (2.*REO+ ] ,+C) ) 

ICODE=0 
GO TO 16 

7 IFIYU ) .i T. (l./(l.+RN + RAR + 2.*Rt.O) ) ) GO TO 16 
Y ( 1 ) =1 ./ ( l . +PN+RAR+2.*RFG) 

16 IF (REQ.EO.0. ) 00 TO 10 
I F ( REO . GT * 1 . ) GO TO IS 

B= (Y (L) *<Y(L> *H12 + R1 1 ) +810) **2-4.*<Y (L) * ( Y ( L > **22 + 42 1 ) + B20) * ( Y (L) 
1* ( Y (L) *602 + 60 1 ) ) 

Z (L) = ( + 1 . *SOPT < B ) - 1 .* (Y <L> * <Y <L) *B12 + B1 1 ) +B10) ) / < ( Y (I.) *< Y ( I . ) * H 2 2 
1 +821 ) +H20) *2. ) 

GO TO 9 

15 B= (Z (L) *(Z(L) *821 +B1 l ) + B01 ) **2-4 . * ( 7 ( L ) * ( / ( L ) *r’22 + H 1 2 ) + H0 2) * (7 (I.) 
1* (7 ( L ) *820 + BIO) ) 

IF (B.LT.O.) GO 70 20 
22 IF (RN.FO. 0 • ) GO TO 21 
IF ( I CODE .EG) • 1 ) GO TO 21 

Y (L) = (-l .*SQRT <R) +1 .* <7 <L> *<7 (L ) *R21 + 811) + B 0 1 ) )/( (Z(l. )*(/(!. )#n?2 + 

1 B 1 2 ) +B02)*(-2. ) ) 

IF ( Y ( L ) • G T . 7 (I ) )G0 TO 2 1 
GO TO 9 

21 Y <L) = < + 1 ,*SQPT (R) +1 .* (7 (L ) * (Z (L) *W?1+H 1 1 ) + 80 l ) ) / ( (/ <1 ) * (/ (I ) *822 + 
1B1 2) + B02) * (-2. ) ) 

ICOOE= 1 
GO TO 9 
10 Z ( L ) =0 . 

9 F (I. > = Y (L> * ( Y (I.) * ( Y (L ) * ( Y <L> *A4u + 7 (L> * (7 (L> *A 32 + A81 ) + A3U ) + 7 <1 ) * ( 

1 7 ( L ) * ( 7 ( L ) * ( 7 ( L ) *A24 + A23 ) + A22 ) + A2 1 ) + A20 ) + 7 ( L ) * ( 7 ( L ) * ( Z ( L ) * ( / ( L ) * 
2A14 + A13)+A1?)+A]1)+A10)+Z(L)*(Z(L)*(7<!.. )*(7(L)*At)4+A03)+An2)+A0i) 

3 + A00 

IF (L.E0.50) GO TO SOO 

OAOY=Y (L) * ( Y (|, ) * ( Y <L ) *4.*A40 + 3«* (7 (|_) * (7 <L) *A32 + A3l ) + A3<' ) ) + 2. * ( 

1Z (L ) * (7 (I ) * (7 (I ) * (7 <L) *A24+A2i> + A22) + A21 ) + A20) ) ♦ / ((_> * (7 (I 1 * (7 (I ) * ( 
2Z (L)*A14+A13) + A 1 2 > + A 1 ) ) +A10 

DAOZ = Y (l ) * (Y (l. ) * ( Y (I. ) * (2.*A32*7 (1 ) + A 11 ) +7 < I. ) * (7 (L) * (7 (I. ) *4.* A 24 + 

1 A23*3. ) +2.*A2?> + A21 ) +7 (L > * </ <L> * (7 <L) *4.*A14+ 3.*A13) +2.*A12) + A1 1 
2) + 7 (L) * (7 (l ) * (7 (L ) *4.*A04+A03* i.)+2.*A02)+A0l 
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DROY=Z ( L ) * ( 7 ( I. )*(?.*R??*Y (I.) + R21 ) +2.*ril2*Y (U + <1 t ) +2 .*m 07*Y (| ) i 
DRDZ=?.*Z <L )*{Y(L)*(Y(I + R20 )+Y(L)*(YiL> + m 1 > + »< 1 

IF (REQ.GT. 1 . ) GO TO 14 

Y (I + ] ) =Y (L ) -F (I ) / (nAr>Y-r)An7*l)nnY/nH;)7 > 

IF ( ( (Y (L+l ) /Y (L) ) .LT. I . 000 1 0) . and. { ( Y ( |_+ l ) /Y ( I ) ) . F f . . murro ) ) GO 10 4 
L = L + 1 
GO TO 16 

14 7 (L + l ) =7 <L) -F (I. > / (DAIl 7-DAMY *1)00/ /nrtOY) 

IF ( ( ICOOF.FO. n ) .AMI). ( (OAD/-0 A' > 7 <H >H0 7/1 ’HI > Y ) . L ] . (! . ) ) GO 10 ,>4 
IF ( ( ICOOF.FO. 1 ) .AMO. ( (DA0/-DAO T'^ijMiiZ/oHOY ) . G T . •) • ) ) GO T (. 24 
IF ( ( (7 (L+l ) /7 (I ) ) .1 T. 1 .0001 0) .AMO. ( ( 7 ( L + 1 > // ( L > > . G T . ,4'Ms-)) ) Gol Ok 
L = L + 1 
GO TO 10 
20 Y (L ) =.Q*Y (L ) 

GO TO 1 6 

26 Z(l ) = .0*7 (I. ) 

GO TO 14 

25 7 (L ) =1 .01 *7 (I ) 

GO TO 14 
20 B=.O01 

PRINT 200 

200 FORMAT (20X27H*** RADICAL I F 4S 1'HAM .) ***> 

GO TO 2? 

27 Z ( L ) =• P*7 ( L > 

GO TO 14 

30 IF (REQ.GT. 1 . ) GO TO 31 
X ( 10) =?.*RFQ/ (C+RFQ) 

X (5) = ( 1 . -KFO ) / (C + RF'O) 

GO TO 11 

31 IF (REQ.GT .20. ) X (7) =1 . 

IF (REQ.GT. 20. ) GO TO l 1 
X ( 10) =?./ (0+?.*RF 0-1 . ) 

X (7) = (REQ-1 . ) »X ( 10) 

X (4) = (4M* (X(l)+?.*X(5)+X(4)+x(D)+y(lO))-(x(4)+x(7)))/^. 

AMK=1 0.** ( (101 14. /T)— 12.044) 

PLIM=. 0 1 / (SQPT ( A M K * X ( 4 ) * X ( 7 ) ** D ) 

IF (P. LI .PL I M) GO TO 11 
PRINT 400 

400 FORMAT ( 10X24H*** AMMONIA HAS FORMFl) **•») 

GO TO 11 
8 X ( 1 ) =Y (L) 

X (4)=Z (L) 

X ( 1 0) =P**?*K1 0*X ( 1 ) *X (4) **2 
X (5) =P*K5*X ( 1 ) **? 

X ( 9 ) =P*KQ*X (4) * X ( l ) 

X (7) =P*K7*X (4) **? 

IF (RN.F Q. 0 . ) GO TO 1 1 

X (2) = (2.- ( 1 . +o) *X ( 1 0 ) -?.*C*X (S) -?.*X ( 7) - ( l . + C) " X ( l ) - ( 1 . +C) *X (P) 

1 -?.*X (4) ) / (C*o»Kft*X ( ) > + 1 . ) 

11 X (8) =P*K8*X ( 1 ) *X (?) 

X (4) = (PN* (X (l)+?.*X(5)+X<4)+X<V)+X(10))-(X(8)+<<?)))/2. 
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X (3) =RAR* ( X ( ] ) +?.*X (S) +X < R) +X ( j ) +X U 0 > ) 

DO 24 J=l«10 

IF ( ( X ( .1) .L T. ft . ) . AMO. (RFO.I.T . 1 . > > GO TO 23 
24 CONTI Ml IE 

PRINT '7 0 0 « ( X ( J ) • J= l » 1 0 ) 

7 00 FORMAT ( RX2HX0. 1 3XPHXM. 1 3X3HXAR * 1 3X2HXH. 13X IHXOP/'-jF lo. m/OX XM2 • 
inX3HX‘-*?« I RXTWXON. 1 3X 3HXOH » 1 RX'tMXHP'VST 10. R / ) 

MM = 0 . 

HM = 0. 

SM = 0. 

DO 47 J=) .10 

IF <X ( J> .FO. 0. ) GO TO 47 
MM=MM+ X (,J)*MR(J) 

HM = HM+ X ( J) <>HR ( .1) 

SM=SM + X ( .1) * ( SP ( ,J) - ALOij ( R <f X < . )) ) ) 

47 CONTINUF 

HM= ( 1 . ORRRR/MM ) *HM 

SM= ( 1 .QRSRR/M k «> *sm 

PRINT R00.MM.HM.SM 

R00 F0RMAT<4XRHMM=T 1 5 . R . ? XRHHM=F i S • R » 7HHTU/L.MM » 2X R>-iSM=t- lS.R.VOMlu/l < ,vl - 
1R//) 

GO TO R00 
END 
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Input 

A single case consists of the determination of the thermodynamic properties of the 
reacting gas mixture, or of one point on a Mollier diagram. It is necessary to input the 
specific heat at constant pressure of each species as a tabular function of temperature, in 
addition to the initial atomic proportions, temperature, and pressure. For the loading 
routine used in the program, any column except the first may be used on the input cards. 
A decimal format is used for the input quantities. A description of the required inputs 
and the names used by the source program are as follows: 


Name 

$NUM1 


HT(1,J) 

GT(1,J) 

CP(I,J) 


Description 

The arbitrary name required by the loading routine to define the input data 
block: specific heat at constant pressure and reference thermodynamic 

properties of each chemical species. 


Reference enthalpy array for each species, nondimensional ; 


Hi 


RT * 


^1 J 

Reference Gibbs free energy array for each species, nondimensional; ’ . 

RT 

Specific heat at constant pressure array for each species, nondimensional; 

Nj.j 


R 

Denotes end of case (column 2). 


$NUM2 

The arbitrary name 
initial proportion ■ 

RN 

r n 

RAR 

r Ar 

REQ 

r eq 

T 

T, °R 

P 

p, atmospheres 


Denotes end of case (column 2). 


The system loading subroutine in the program (NAMELIST) is quite flexible in that 
the order of the input cards is unimportant and successive cases can be run by repeating 
the identification and $NUM cards followed by only the changed parameters and a $ card. 
An example of a set of input cards is given by the following listing of the inputs necessary 
to compute the thermodynamic properties of the stoichiometric hydrogen-air mixture at 
a temperature of 5000° R (2778° K) and a pressure of 1 atm. 


19 



APPENDIX A 


$NUM 1 HT ( 1 * ) ) = 269. 79, HT ( 1 , 2 ) =5 1 2 . Oft * HT ( 1 * 3) =2 . 5003 , HT ( 1 . 4 ) =236 . 28 , 

HT ( 1 ,5) =3.489R,HT ( 1 . 6 ) = 3 . 4853 , Hy ( 1 , 7) =3.7283*HT ( 1 , 8 ) = 1 0 0 . 99 , H f (1 , 9 ) =45.3 22. 

HT ( 1 , 10) =-254.66.GT (1*1) =253.1 2, GT ( 1 ,2) =496 . 1 0 , GT ( I . 3 ) =- 1 3 . 666 . GT ( 1 ,h) = 224.98 , 
GT ( 1 ,5) =-17.746,GT ( 1 ♦ 6) = -l 6 . 1 20 * GT ( 1 * 7 ) =-8 . 8676 * GT ( 1 *8) =79.2 73, GT ( 1 ,9) = 27. n 1 I . 
GT<1.10>=-273.44, 

CP ( 1 , 1 ) =2 • 8 374 • 

?. 7 7 IS *2. 71 0 7. 2.652;?, 2. 61 20,2.5884,2.5690,2.5568,2.5471 ,2.539 1. 
2.5336,2.5288,2.5249,2.5218,2.5193,2.5174.2.5135.2.3141 .2. 51 28. 2. 51 1<S, 2. 51 06. 
2.5098,2.50 9 0.2.5081 , 2 . 50 75 , 2 . 50 7 0 , 2 . 50 68 , 2 . 5064 , 2 . 5057 . 2 . 5034 . 2 . 3053 . 2 . 505 3 . 

2. 5050, 2. 5049, 2. 5050,2.5049, 2. 5050, 2. 5053 .2.3036, 2. 3059, 2, 3 051,2.5064,2. 5 <>/()• 
2.5080 ,2.5084,2.5092,2.5103,2.51 1 4 , 2 . 51 25 . 2 . 5 1 38 , 2 . 3 1 52 . 2 . 5 1 66 , 2 . 5 1 80 , 2 . 3 ] 9 7 . 
2.5214,2.5234,2.5254,2.5277,2.5299,2.5322,2.5344,2.3367,2.3391 ,2.541 7.2.5443. 
2.5473,2.550 1 .2.5530.2.5562,2.5591 , 2 . 56 1 8 , 2 . 5632 . 2 . 3685 , 2 . 57 1 5 . 2 . 574 7 . 2 . 378 1 . 
2.5814, 2.5848, 2. 5882, 2. 5913.2.5949. 2.3982, 2.60 16,2.6049,2.6083,2.61 16.2.6130 . 
2.6183,2.621 7,2.6251 .2.6284,2.631 8,2.6351 ,2.6383,2.6418,2.6432,2.6484,2.6314, 
2.6544,2.6576,2.661 1 , 2 . 6642 . 2 . 6668 , 2 . 6693 , 2 . 6 729 , 2 . 6758 , 

CP ( 1 ,2) =31*2.4999, 

2.5 0 00, 2.50 04, 2. 50 03, 2. 50 04, 2. 500 7, 2. 5010, 2. 50] 3, 2. 50 IS. 2. 50! 7. 
2, 5023, 2.5033. 2 .5041 ,2.5050 ,2.5060,2.5072,2.5086,2.5101 ,2.51 1 7,2.5137,2.5161 , 
2.5185,2.5211 ,2.5244,2.5278,2.5313 .2. 5852. 2. 5394 .2. 3439,2. 3487, 2. 5539, 2. 5593* 
2.5650,2.5711 . 2 . 5775 , 2 . 584? , 2 . 59 1 1 ,2.5986,2.6065,2.6147,2.6232,2.631 7,2.6406, 
2.6501 ,2.6599,2.6700,2.6802,2.6905,2.7013,2.7123.2. 723 7 , 2 . 7352 , 2 . 7469 . ? . 7589 . 
2.7712,2.7837.2.7963,2.8091 , 2 . 8220 , 2 . 835 1 , 2 . 8483 , 2 . 86 1 3 , 2 . 8 747 , 2 . 8885 , 2 . 9022 . 

2. 9 1 58, 2. 9297, 2. 9438, 2. 9577, 2. 97 17. 8. 9R5 7. 2. 9997. 3. 01 36, 3.0276,3.0416,3.0556. 

CP ( 1 ,3) =106*2.4999, 

CP ( 1 ,4) =106*7.4999, 

CP ( 1,5) =3.5001 , 

3.4995,3.5077,3.5223,3.5568,3.6099.3.6698, 

3.7393,3.8077,3.8721 ,3.9329,3.9883,4.0386,4.0841 ,4.1246,4. 1616,4.1947, 
4.2244,4.2518, 

4.2767.4.2999,4.3213,4.34] 7,4.3610,4.3795.4.3970,4.4146,4.4315, 
4.4477,4.4639,4.48 0 1 ,4.4962,4.5] 21 , 4 . 5278 , 4 . 54 34 , 4 . 3388 , 4 . 5 742 . 4 . 5896 , 4 . 6049 , 
4.6204,4.6336.4.650 6,4.6656,4.68 0 3.4.605 1 . 4 . 7099 . 4 . 724 1 . 9 . 7383 . 4 . 7523 • 4 . 766 0. 
4.7793,4.7927,4.8061 ,4.8188.4.8312,4.0436,4.855 7 . 4 . 3674 . 4 . 8 790 . 4 . 8904 , 4 . 90 1 3 . 

4. 9 123, 4. 92?7.4.9329,4. 9<+ 31 , 4 . 93?9 , 9- . 0624 , 4 . 97 1 6 , 4 . 9803 , 4 . 9892 . 4 . 99 78 . 3 . 0 03 7 . 
5.0135,5.021 0,5. 0285, 5. 0358, 5. 0427. 3. 0494, 3. 0559, 5. 0622. 3. 0681 ,5.0740.5.0804. 
5.0861 ,5. 0909,5.096 0 ,5. 1 010.3. 1 054.3, 1 ] 06.3. 1 151 ,3. 1 193,5. 1 237,3. 1 278.6.1316. 

5. 1 355,3.1 395,5. 1482,5. 1468,3. 1499,3. ] 529,3.1560,3.1391,5. 1622.5.1630,5. 16 74. 
5.1699, 

CPU ,6>=3.5003, 

3.50 05,3.5011,3.5012,3.5064,3.5130, 3.3 323,3.53 7 7, 3 . 39 ] I , 3 . 629 o . 

3. 6727, 3. 7169. 3. 76?7,3.8072, 3. -851 ] , i. 8928 . 3 . 932b , 3 . 9698 ,4. OI)<tH ,4 . 037 3 . 8 . 06 78 . 
4.0964.4.1227,4. 1468,4. 1695,4. 1907.4.21 02 . 4 . 2285 . 4 . 2434 , 4 . 26 1 3 , 4 . 276 l . '4 . 29 0 0 . 
4.3029,4.3151,4. 1266,4.3374,4. 3477.4.3573,4.36^3.4. 8 796 . 4 . 382 7 , 4 . 39(>5 . 4 . 39 7H. 
4.4046,4.41 1 0 ,4.41 74,4. 4 236,4. 4295.“+. 4351 . 4 . 44 04 , 4 . 4434 . 4 . 43 0 0 . 4 . 4344 .4 . 4590 • 
4.4633,4.46 71,4.471 1 ,4. 4754. 4. 4793, 4. 4829*4.48 6 3. 4. 4896.4.993 0.4. 49 6 3, 9.9498, 

4. 5029, 4. 5035.4. 5085, 4. 511 7. 4. 5143.4. 5 173. 4. 320 3«4. 32 3 3, 4. 3234, a. c -27 8, +.5 30 8, 
4.5329,4.5354,4.5377,4.5 399.4.5421 , 4 . 5444 , 4 . 54 7 1 ,4.8493.4.331 2,8.3332.9.3341. 

4. 5572, 4. 5595, 4. 5633, 4. 5714, 4. 5737. h. 5761 , 8 . 3 782 . 4 . Sb0 6 . 4 . -,832 , 4 . 5 83 7 . 4 . 38 h2 . 
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4.5907,4.5933.4.5962,4.5988,4.6013.9.6039,4.6068.4.6096. 

CP(1,7)=2.7955, 

3.133) .3. 3547. 3. 4404. 3. 4995. 3. 6064. 3. SI 73. 

3. 51 89, 3. 5234. 3. 5?77. 3. 5346. 3. 5442. 3. 5561 .3.5709, 3.5885, 3 . 609 1 , 3 . 6 1?6 . 
3.6577,3.6837,3.7107,3.7384,3.7670.3.796) ,3,8260, 3, h 655,3.8847 ,3.913 / « 3.9'*??. 
3. 9703, 3. 9978, 4. 0247,4. 0506,4. 0756, 4. 1000, 9. 1238, 4. 1 466 .4. 1 690 ,4 . 1 9 1 0 ,^. v» l 22 . 

4.2324.4.2525.4.2724.4.2915.4.3100.4.3280.4.3446.4.8628.4. -1795,4.3966,4.4) ) 6. 
4.4271 ,4.4426,4.4579,4.4728,4.48 78.4.5021 , 4 . 5 1 64 . 4 . S3 0 3 , 4 . 8448 , 4 . SSH .1 , /, . 87 ? 1 . 
4.5857,4.5989,4.6120,4.6252,4.6 382.4.651 1 , 4 . 6638 , 4 . 6763 . 4 . 688 7 . 4 . 70 09 . 4 . 7 1 6. . 

4. 7255. 4. 7372. 4. 7489. 4. 7607. 4. 7724. 4. 7842. 4. 7987. 4. HO 71 ,4.8184,4.8298,4.840 1. 
4. 8512. 4. 8621, 4. 8730*4. 8840.4. 8947. 4. 9052. 4. 9156, 4. 9259, 4 . 9 3b 3. 4. Q466.4.98 70. 

4.9672.4.9772.4.9871 , 4 . 9968 . 5 . 0067 , 5. 0 1 64 , 5 . 0288 , 8. 0354 , 5 . 0453 , 5 . 084 3 , 8 . 04 39 . 
5.0731 , 

CP ( 1 ,8) =3.8805. 

3.7195,3.6342.3.6006,3.5865,3.5991 ,3.4252.3.6679,3.7173,3. 7482. 

3.8216.3.8728.3.9227.3.9691 .4.0124.4.0522.4.0888.4. 1217.4.1823.4. ) HO 1 ,4.9087. 
4.2291 ,4.2507,4.27 0 5,4.2886.4. 3 049.4.3222,4.33 75 . 4 . 350 3 , 4 . 3624 , 4 . 3746 , 4 . 3°o h . 
4.3932,4.4022,4.411 1 .4.4197,4.42/8,4.4354.4 .4426 ,4 .449? . '♦.♦♦558 , 4 . 4622 . 4 .4483 , 
4.4740 ,4 .4 794. 4.4R46, 4 . 4896,4.4946,4.4997,4.8045.4.8069,4.51 3 1.0.8171. 4. 52 15, 
4.5256,4.5295,4.5334.4.5374.4.541 1 ,4.544 7, 4. 54 84, 4. 552 0.4. 5550,4.55 87, 4. 5621 . 
4.5654,4.5688,4.5721 ,4.5758,4.5787,4.8817,4.585 1 , 4 . 5h84 . 4 . 59 1 9 , 4 . 5945 , 4 . 59 76 • 
4. 6006, 4. 60 3 7, 4. 6068, 4. 60 99, 4. 6 126. 4. 6 158, 4. 61 86 «4. 6 21 7, 4. 0248,4.627 7. 4. 63 04. 
4.6334,4.6365.4.6393,4,64?] ,4, 6449, ‘+.6477, 4. 65 08, 4. 68 34, 4. 6846 , 4 . 6898, 4. 66 ??. 
4.6647*4.6673.4.6701 « 4. 6729, 4. 6757, 4. 6784, 4. 6808,4.6835. 

CP ( 1,91=3.704?, 

3.6653,3.6321 ,3.6051 ,3.8823.3.5640,3.8426, 

3.5476.3.5460.3.550 7,3.8591 . 1 . 57 1 9 , 3 . 5895 . 3 . 6 1 0 3 . 3 . 6 346 . 3 . 66 0 9 , 3 . 6895 . 
3. 7 194, 3.75 00.3. 781 5,3.81 2 7. 3.84 38,3.8743. 3.9 045,3. 9338, 3.9622, 3. 99 0 0, 0 . 0 168. 
4.0428,4.067 7.4.09 1 7, 4. 1150, 4.1 3 75«‘+. 1590, 4. 1796, 4. I 995.4.2186,4.2368,4.2546 . 
4.2720,4.2886,4. 3043,4.3197.4.3346.4. 3490.4.3631 ,4. 1769,4.3902,4.4030,4.4186. 
4.4278,4.4397.4.451 4,4.4629.4.4742.4.4851 ,4.4959,4.8068,4.81 7 0,4.5273,4.8378, 
4.5475,4.5573.4.5670,4.5769,4.5868,4.8960,4.6055,4.6150.4.6246.4.6340,4.6436. 
4. 65 30, 4. 66?]. 4. 67 15. 4. 68 10, 4. 6905. 4. 700 1,4. 7096, 4. 7 191. 4 . 7286.4.7381,4.7476. 
4.7571 ,4.7665,4.7762,4.7861 .4.7957,4.8051 .4.7960,4. 7943,4.8029,4.8099,4.8165. 
4.8232,4.830 0 ,4. 8363. 4. 84? 3. 4. 848 3. 4-. 8541 ,4.8597,4.8661 ,4.8706,4.8786.4.8802. 
4.8847. 

CP<1, 101=4.0048, 

4.0059,4.0143,4.0279,4.0598,4.1086.4.) 662, 
4.2345,4.3064,4.3809,4.4583,4.5377.4.6192,4. 7621 ,4.7867,4.8714,4.9666. 
5.0412,5.1246,5.2071 ,5.2873,5.3654,6.4409,5.51 39 , 8 . 5844 . 5 . 6528 , 8 . 7 1 76 , 8 . 780 0 . 
5.8394,5.8964.5.950 9,6.0 031 ,6.0 529.6. 1005. 6. 1456, 6. 1889. 6. 23 0 3. 6. 269 7. 6. 10 75. 
6.3435,6.3780,6.41 ) 0 , 6 . 4424 . 6 . 4722 , b . 5 0 1 2 . 6 . 5290 . 6 . 6554 , 6 . 58 1 0 , 6 . 6056 , 6 . 629 1 . 
6.6515,6.6732,6.694) ,6.7143.6.7337,6.7522,6.7703,6. 788 0 . b . 8048 , 6 . 8208 . 6 . 8 365 . 
6. 8516.6. 8663, 6. 8807, 6.8948, 6. 9083. 6. 9?1 0.6. 9337.6. 9464, 6. 958 0,6. 9689. 6. 9 h 06. 
6.99)7,7.0019.7.0] 23 , 7 . 0228 . 7 . 0326 . 7. 04 15, 7. 0509, 7. 0605. 7. 068 7. 7. 077 3. 7. 0866. 
7.0949,7.1022,7.1100,7.1 183,7.1257, 7 . 1 324 , 7 . 1 377 , 7 . 1 43 1 , 7.1484,7.1836.7. 1586. 
7. 1637,7. 1691 .7. 1743,7. 1792,7. 1 846, 7.1898,7.1948*7. 1999, 7 . 2049 . 7 . 20 99 . 7 . 2 1 50 . 
7.22005 

SNUM2 RN=3. 727445, RAR=. 02296952, REO= 1. 0 ♦ T=50 00.0, P=. 00 18a 
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APPENDIX A 


It should be noted that the input data in block $NUM1 were obtained from refer- 
ence 2 and need not be changed from case to case unless fewer points are desired. The 
present coverage includes temperatures from 200° R to 10 000° R (111 0 K to 5556° K) 
in 100° R (56° K) increments. 


Output 

The output names used in the source program and the corresponding quantities from 
the main body of the report are as follows: 

Program name Quantity 


xo 

x o 

XN 

X N 

XAR 

X Ar 

XH 

X H 

X02 

x o 2 

XN2 

x n 2 

XH2 

x h 2 

XON 

X ON 

XOH 

X OH 

XH20 

x h 2 o 

MM 

M 

HM 

H 

SM 

S 
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An output listing for the example input is as follows: 

RN=3. 72744 500E+00 RAR=2.29695200E-02 REQ=1.00000000E+00 



T=5.00000000E+03R P=1.00000000E+00ATM 


xo 

5.70327244E-03 

XN 

2.31249727E-06 

XAR 

7.60417436E-03 

XH 

1.61297775E-02 

X02 

1.34353754E-02 

XN2 

6.12842267E-01 

XH2 

4.58018477E-02 

XON 

8.30244316E-03 

XOH 

2.59805399E-02 

XH20 

2.64197990E-01 


MM=2.35537499E+01 HM=7.96227455E+02BTU/LBM SM=2.76587239E-tOOBTU/LBM-R 
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APPENDIX B 


FORMULATION OF INITIAL VALUES FOR THE NEWTON ITERATION SCHEME 


The purpose of this appendix is to present the simplified gas models which were 
used to estimate the desired solution over a wide range of variables. To make certain 
that the simplified gas model gave results nearly equal to the entire solution, different gas 
models were chosen in separate regions of equivalence ratio. 

For r eq < 1: 

Gas model chosen 


H2O, O2, O, N2, Ar 
Conservation of mass equations 

x H 2 0 + X 0 2 + X 0 + %2 + X Ar = 1 


2X N 2 = r N( X H 2 0 + 2X 0 2 + X o) 
X Ar = r Ar( X H 2 0 + 2X 0 2 + X o) 
2X H 2 0 = 2r eq( X H 2 Q + 2X Q 2 + X o) 


Equilibrium expression 
X 0 2 = P K 5 X 0 
Solving for Xq gives 


v _ y(l + C) 2 +16(1 - r eq )(C + r eq )pK 5 - (1 + C) 

X ° = 4(C + r^pKj 

The presence of an extraneous solution to the entire problem made it necessary to 
establish a correction factor which would insure convergence to the proper solution. The 
correction factor was determined by comparing beforehand the result of the simplified 
gas model with the entire solution over the range of variables. The correction factor 
chosen for r eq < 1 is 


x o(D 

x o 


= 1.1 - (0.5 X 10“ 4 )t 
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The initial estimate for the iteration becomes 


Xq(D = 


|/(1 + C)^ + 16^1 - r e q^C + r eq)P-^5 - ^ + 


4(C + r eq )pK 5 


x 





X 10" 



For r eq = 1: 

Gas model chosen 

H 2 0, O, H, N 2 , At 
C onservation of mass equations 

x H 2 0 + x O + X H + X N 2 + x Ar = 1 

2x N 2 = r N( x O + X H2 o) 

x Ar = r Ar( x O + x H 2 d) 

X H + 2X H 2 0 = 2r eq^ x O + X H 2 o) 

Equilibrium expression 

x h 2 o = P 2 K 10 Xcp 4 
Solving for Xq yields 

X 0 * (2(1 + C)p2K 10 ] 1/3 
The correction factor chosen for r eq = 1 is 


X °^-- = 0.34 - (0.2 x 10 " 4 )t 
x O 


The initial estimate for the iteration becomes 


Xq(1) = 
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For r eq > 1: 

Gas model chosen 

H 2 , H, O, N 2 , Ar 

Conservation of mass equations 

x H 2 + X H + x O + x N 2 + XAr = 1 

2x N 2 = r N x O 

x Ar = r Ar x O 


2Xh 2 + X H = 2r eq X 0 
Equilibrium expression 
x H 2 = P k 7 x H 

Solving for Xjj and determining that no correction factor was required for this case 
yields the following initial estimate: 

_ \Jfi r eq + 1 + + 32r eq pK7(2r e q + 1 + C) - (4r e q + 1 + c) 

H ~ 4(2r eq + 1 + C)pK ? 


26 



REFERENCES 


1. Zeleznik, Frank J.; and Gordon, Sanford: An Analytical Investigation of Three General 

Methods of Calculating Chemical -Equilibrium Compositions. NASA TN D-473, 1960. 

2. Browne, W. G. ; and Warlick, D. L. : Properties of Combustion Gases — System: 

H 2 -Air. R62FPD-366, Advan. Engine Technol. Dep., Gen. Elec. Co., Nov. 1962. 
(Reprinted Sept. 1964.) 

3. Erickson, Wayne D. ; Kemper, Jane T.; and Allison, Dennis O.: A Method for Computing 

Chemical -Equilibrium Compositions of Reacting-Gas Mixtures by Reduction to a 
Single Iteration Equation. NASA TN D-3488, 1966. 

4. Mascitti, Vincent R.: A Simplified Equilibrium Hydrocarbon-Air Combustion Gas 

Model for Use in Air-Breathing Engine Cycle Computer Programs. NASA 
TND-4747, 1968. 

5. Franciscus, Leo C. ; and Lezberg, Erwin A. : Effects of Exhaust Nozzle Recombination 

on Hypersonic Ramjet Performance: H. Analytical Investigation. AIAA J., vol. 1, 
no. 9, Sept. 1963, pp. 2077-2083. 

6. McCracken, Daniel D.: A Guide to FORTRAN Programming. John Wiley & Sons, Inc., 

c. 1961. 


27 



to 

CO 


TABLE L- CHEMICAL COMPOSITION AND THERMODYNAMIC PROPERTIES OF AIR (r eq = o) 
[First value of paired numbers, present method; second value, ref. z\ 



T = 

2000° R (1111° K) 

T = 

6000° R (3333° K) 

T = 10 000° R (5556° K) 


p = 0.001 atm 

p s= 1.0 atm 

p = 100 atm 

p = 0.001 atm 

p = 1.0 atm 

p « 100 atm 

p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

x o 

0 

0 

0 

0.34007 

0.10807 

1.2953 x lO' 2 

0.22018 

0.32149 

0.22343 


4.4792 x 10-8 

1.5018 X 10-9 

0 

.34011 

.10836 

1.2995 X 10-2 

.22014 

.32154 

.22437 

x N 

0 

0 

0 

2.436 X lO' 3 

8.1153 x lO' 5 

8.2882 X lO' 6 

0.72846 

8.1383 x lO" 2 

8.5891 x lO’ 3 

0 

0 

0 

2.439 X lO' 3 

8.1193 x lO’ 5 

8.2971 X lO' 6 

.72869 

8.1557 X 10' 2 

8.6055 X 10' 3 

^Ar 

9.624 X XO’ 3 

9,624 x lO' 3 

9.624 X 10-3 

7.9758 X 10-3 

9.1036 X 10“ 3 

9.5616 x 10-3 

0.50591 X lO' 2 

7.6853 X lO' 3 

8.5075 X lO* 3 


9.624 x lO- 3 

9.624 X lO' 3 

9.624 x lO' 3 

7.9757 X 10’ 3 

9.1022 x 10" 3 

9.5614 X 10" 3 

.50583 x lO' 2 

7.6843 X lO" 3 

8.5029 X lO -3 

X H 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

X Op 

0.20949 

0.20949 

0.20949 

1.1752 X 10 -3 

0.1186 

1.7049 X lO’ 1 

0.291429 X lO -6 

6.2133 X 10' 4 

3.0010 X 10' 2 

2 

.20944 

.20944 

.20944 

1.1684 X 10‘ 3 

.1186 

1.7056 X lO -2 

.28385 X lO -6 

6.0554 X 10' 4 

2.9485 x 10 -2 

% 2 

0.78088 

0.78088 

0.78088 

0.64353 

0.71315 

0.74464 

4.6227 X lO" 2 

0.57696 

0.64253 


.78083 

.78083 

.78083 

.64353 

.71314 

.74472 

.6047 X lO -2 

.57682 

.64219 

X H 2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X ON 

0 

0 

0 

4.8120 X 10‘ 3 

5.0931 X lO" 2 

6.2348 X lO -2 

0.72654 X lO' 4 

1.1851 X lO’ 2 

8.6929 X 10‘ 2 

1.0302 X 10-4 

1.0302 X 10-4 

1.0302 X lO' 4 

4.7821 X lO' 3 

5.0719 X 10' 2 

6.2156 X lO -2 

.72157 X 10' 4 

1.1796 X lO -2 

8.6848 x 10 -2 

X OH 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

x h 2 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

X NH 

l 

0 

0 

0 

0 

0 

0 

0 

0 

0 

x NH 3 

0 

0 

0 1 

0 

0 

0 

0 

0 

0 

M 

28.965 

28.965 

28.965 

24.005 

27.399 

28.777 

15.227 

23.130 

25.605 


28.965 

28.965 

28.965 

24.004 

27.395 

28.777 

15.224 

23.128 

25.591 

H, Sis. 

504.7 

504.7 

504.7 

3296.1 

2210.2 

1835.5 

14 687.4 

5275.8 

4135.6 

’ lbm 

504.8 

504.8 

504.8 

3296.3 

2211.2 

1835.5 

14 692.4 

5278.6 

4141.2 

H, JL 
’ kg 

1.1731 x 10® 

1.1731 X 10® 

1.1731X108 

7.6615 X 10® 

5.1374 X 10® 

4.2664 X lof 

34.139 X 10® 

12.263 X 10§ 
12.270 X 10 6 

9.6128 X 10® 

1.1734 X 10® 

1.1734 X 10 6 

1.1734 X 10 6 

7.6619 x 106 

5.1397 x 10 6 

4.2664 X 10 6 

34.151 X 10 6 

9.6258 X 10® 

s Btu 

2.4444 

1.9708 

1.6551 

3.1206 

2.3941 

2.0084 

4.4489 

2.7886 

2.2974 

' lbm-°R 

2.4445 

1.9709 

1.6552 

3.1208 

2.3944 

2.0084 

4.4497 

2.7890 

2.2979 


10 227 

8245.8 

6924.9 i 

13 057 

10 017 

8403.1 

18 614 

11 668 

9612.3 

' kg-°K 

10 228 

8246.2 

6925.4 , 

13 057 

10 018 

8403.1 

18 618 

11 669 

9614.4 







TABLE n.- CHEMICAL COMPOSITION AND THERMODYNAMIC PROPERTIES OF HYDROGEN-AIR COMBUSTION FOR r art = 1 
[First value of paired numbers, present method; second value, ref. 2j 


T = 2000° R (1111° K) T = 6000° R (3333° K) T = 10 000° R (5556° K) 



p = 0.001 atm 

p - 1.0 atm 

p = 100 atm 

p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

Xo 

0 

0 

0 

0.20125 

4.8479 x 10" 2 

2.8052 x lO' 3 

0.15156 

0.19233 

0.12007 


0 

0 

0 

.20126 

4.8651 X lO" 2 

2.8201 x 10-3 

.15156 

.19238 

.12064 

X N 

0 

0 

0 

1.8725 x lO" 3 

6.9610 X lO' 8 

7.5493 x 10" 6 

0.51841 

6.2263 X lO- 2 

6.9426 X lO' 3 


0 

0 

0 

1.8750 x 10-3 

6.9635 x lO -5 

7.5595 x HT 8 

.51843 

6.239 x 10“ 2 

6.9498 X 10" 3 

*Ar 

7.9570 X 10- 3 

7.9570 x 10“3 

7.9570 x 10' 3 

4.7123 x 10-3 

6.5887 x lO" 3 

7.6958 x lO" 3 

3.4821 x lO- 3 

4.5792 X lO* 3 

5.4487 x lO" 3 


7.9570 x 10“ 3 

7.9570 x lO" 3 

7.9570 x 10~ 3 

4.7123 x 10" 3 

6.5859 x 10 “ 3 

7.6948 x 10" 3 

3.4821 x lO" 3 

4.5793 x 10' 3 

5.4478 x lO" 3 

X H 

0 

0 

0 

0.40729 

0.11748 

7.8295 x lO" 3 

0.30319 

0.39507 

0.27773 


3.1673 x 10~ 9 

0 

0 

,40729 

.11768 

7.8495 x 10" 3 

.30319 

.39487 

.27709 

x o 2 

0 

0 

0 

4.1154 x lO" 4 

2.3883 x lO' 2 

7.9961 x lO" 3 

1.3808 x 10‘ 7 

2.2237 x lO" 4 

8.6659 x lO' 3 


4.2358 x lO’ 6 

3.9609 x lO' 7 

7.6624 x 10" 8 

4.0914 x 10- 4 

2.3907 x lO" 2 

8.0326 x lO" 3 

1.3457 x 10- 7 

2.1677 X 10~ 4 

8.5250 X lO* 3 

x n 2 

0.64562 

0.64563 

0.64563 

0.38032 

0.52477 

0.61827 

2.3315 x lO' 2 

0.33771 

0.41975 


.64562 

.64563 

.64563 

.38032 

.52456 

.61819 

2.3307 X lO" 2 

.33756 

.41886 

x h 2 

0 

0 

0 

1.0593 x 10“ 3 

8.8132 x lO' 2 

3.9147 x lO" 2 

7.3389 x lO" 7 

1.2460 x lO- 3 

6.1580 x lO -2 


9.1818 X lO" 6 

9.4953 x lO' 7 

2.1589 x 10“ 7 

1.0573 x lO” 3 

8.8271 x 10" 2 

3.9272 x lO” 2 

7.3418 x lO' 7 

1.2453 X 10" 3 

6.1322 x lO" 2 

x ON 

0 

0 

0 

2.1885 X lO" 3 

1.9599 x 10~ 2 

1.2299 X 10‘ 2 

3.5591 x 10' 5 

5.424 x 10‘ 3 

3.7759 x 10“ 2 


4.2126 x lO' 7 

1.2882 x lO’ 7 

5.6659 x 10" 8 

2.1755 X 10-3 

1.9530 x lO" 2 

1.2290 x lO' 2 

3.5343 x lO" 5 

5.3989 X lO" 3 

3.7715 X 10“ 2 

XOH 

0 

0 

0 

8.9308 x lO" 4 

6.205 x lO- 2 

2.3931 x 10-2 

6.9703 X lO' 7 

1.1525 X lO" 3 

5.0582 X lO- 2 


5.8048 x 10- 7 

5.7083 x 10-8 

1.1971 x 10" 8 

8.9124 X lO' 4 

6.224 x lO" 2 

2.4067 x lO" 2 

6.9554 x 10' 7 

1.1498 X 10 -3 

5.0599 x lO” 2 

x h 2 o 

0.34642 

0.34642 

0.34642 

5.4362 x 10-8 

0.10895 

0.28002 

1.7262 x lO" 12 

3.7192 x 10’ 6 

1.1474 X 10' 2 


.34641 

.34642 

.34642 

5.3754 X 10‘ 6 

.10848 

.27976 

0 

3.6164 x 10‘ 6 

1.1168 x lO -2 

X*TTT 










NH 

0 

0 

0 

1.9938 X 10 -6 

2.1395 x 10 -5 

1.5492 x HT 5 

1.3726 x 10- 6 

2.1514 x 10" 4 

1.6817 x 10' 3 

x nh 3 


— 

— 

--- 



--- 

— 

— 


0 

0 

0 

0 

9.5618 x lO' 8 

3.0803 x 10" 6 

0 

0 

0 

M 

24.647 

24.647 

24.647 

14.596 

20.408 

23.837 

10.785 

14.184 

16.877 


24.647 

24.647 

24.647 

14.596 

20.400 

23.834 

10.786 

14.184 

16.874 

H Btu 

-814.9 

-814.9 

-814.9 

6637.4 

2650.5 

1115.9 

18 610.5 

9314.3 

6651.4 

’ Ibm 

-814.8 

-814.9 

-814.9 

6637.6 

2655.3 

1117.3 

18 611.1 

9317.9 

6667.9 


-1.8942 x 106 

-1.8942 X log 

-1.8942 x 10 8 

15.428 X 10 8 

6.1608 x 10 8 

2.5938 x 10 8 

43.258 x 10 8 

21.650 X 10 8 

15.461 X 10 8 

H ’ kg 

-1.8939 x 10 8 

-1.8942 X 10 6 

-1.8942 x 10 6 

15.428 x 10 8 

6.1720 x 10 6 

2.5971 x 10 6 

43.260 X 10 8 

21.659 X 10 8 

15.499 X 10 6 

s Btu 

2.8702 

2.3136 

1.9425 

4.6147 

3.0999 

2.4381 

6.0321 

3.9998 

3.1268 

5 lbm-°R 

2.8701 

2.3135 

1.9424 

4.6147 

3.1008 

2.4383 

6.0325 

4.0004 

3.1287 


12 009 

9680.1 

8127.4 

19 308 

12 970 

10 201 

25 238 

16 735 

13 083 

S ’ kg^K 

12 008 

9679.7 

8127.0 

19 308 

12 974 

10 202 

25 240 

16 738 

13 090 



NASA-Langley, 1969 33 L-6627 


Co TABLE m.- CHEMICAL COMPOSITION AND THERMODYNAMIC PROPERTIES OF HYDROGEN-AIR COMBUSTION FOR r eq = 2.5 


First value of paired numbers, present method; second value, ref. fj 



T = 

2000° R (llll c 

>K) 

T = 

6000° R (3333° K) 

T = 

10 000° R (5556 

i° K) 


p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

1 p = 0.001 atm 

p = 1.0 atm 

p = 100 atm 

X Q 

0 

0 

0 

0.11621 

1.9146 x 10“2 

2.577 x lO" 4 

0.10372 

0.12046 

3.1192 x lO" 2 


0 

0 

0 

.12498 

1.9260 X 10- 2 

2.6015 x 10- 4 

.10372 

.12049 

' 7.1178 x lO* 2 

X N 

0 

0 

0 

1.2330 x 10~ 3 

5.7712 X 10“ 5 

6.2398 x 10' 6 

0.36368 

L . .8672 x 10-2 

6.4016 X lO’ 2 


0 

0 

0 

1.4765 X 10-3 

5.7675 x 10~ 5 

6.2149 X 10-3 

.36372 

4.8747 x 10-2 

5.6495 x 10" 3 

X Ar 

5.2362 X 10" 3 

5.2362 X 10~ 3 

5.2362 X KT 3 

2.7139 x 10- 3 

4.4766 X lO" 3 

5.1558 x 10“ 3 

2.3829 x lO" 3 

2.8578 X lO' 3 

3.4546 X lO" 3 


5.2362 X 10-3 

5.2364 X 10“3 

5.2566 x 10~ 3 

2.9222 x 10-3 

4.4749 X 10" 3 

5.1559 x 10' 3 

2.3828 x lO' 3 

2.8580 x 10" 3 

3.5702 x 10“ 3 

X H 

0 

0 

0 

0.65611 

0.19685 

2.2727 x lO" 2 

0.51869 

0.61491 

0.43121 


6.1123 x 10~ 7 

1.9327 X 10~ 8 

1.9201 x 10" 9 

.63017 

.19705 

2.2746 x lO" 2 

.51869 

,61470 

.41915 

x o. 

0 

0 

0 

1.3723 X lO" 4 

3.7249 x 10- 3 

6.7485 X 10" 5 

6.4672 x 10" 8 

8.722 x lO* 5 

5.8487 x lO' 4 


0 

0 

0 

1.5776 X 10' 4 

3.7467 X 10“ 3 

6.8360 X 10" 5 

6.3018 x lO' 8 

8.5033 X lO” 8 

2.9673 x 10-3 

X No 

i 0.42486 

0.42486 

0.42486 

0.21917 

0.35999 

0.41786 

1.1492 x lO” 2 

0.20621 

0.20307 

iN 2 

.42486 

.42486 

.42457 

.23583 

.35984 

.41783 

1.1472 X lO" 2 

.20607 

.27676 

X H ? 

0.34194 

0.34190 

0.34194 

2.7491 X 10" 3 

0.24745 

0.32985 

2.1479 x lO" 8 

3.0186 X 10-3 

0.14844 

z 

.34194 

.34190 

.33744 

2.5311 x 10" 3 

.24748 

.32977 

2.1487 x lO"® 

3.0179 X lO' 3 

.14032 

X 0N 

0 

0 

0 

8.3215 x 10' 4 

6.4171 x 10' 3 

9.3387 X lO' 4 

1.7087 x lO' 5 

2.6557 x 10“ 3 

9.0449 x lO' 2 


0 

0 

0 

1.0638 x 10" 3 

6.4037 x 10' 3 

9.3207 x 10" 4 

1.6970 x lO' 5 

2.6420 X lO' 3 

1.8088 x lO- 2 

X OH 

0 

0 

0 

8.3077 X 10" 4 

4.1065 x 10-2 

6.3816 x 10" 3 

8.1607 x 10" 7 

1.1235 X lO' 3 

2.0402 X 10-2 


1.9795 x 10-9 

0 

0 

8.5629 X 10” 4 

4.1262 X lO’ 2 

6.4338 x lO" 3 

8.1431 X 10- 7 

1.1211 X lO' 3 

4.5157 x lO" 2 

X H«0 

0.22796 

0.22796 

0.22796 

8.1463 X lO' 6 

0.12081 

0.21675 

3.457 x 10- 12 

5.6430 X lO’ 6 

7.1857 x 10-2 


.22796 

.22797 

.22885 

7.9908 x lO' 6 

.12041 

.21671 

0 

5.4891 x lO" 6 

1.5077 x lO' 2 

X NH 

— - 

--- 

- — 

— 

— 

— 

... 

— 

--- 


0 

0 

0 

2.4292 x 10-3 

2.9671 x 10“ 5 

3,6907 x 10~ 5 

1.6475 X lO' 6 

2.6168 x lO' 4 

2.0679 x 10~ 3 











x nh 3 

3.9734 x 10~ 8 

3.9726 x 10-5 

3.8938 x 10-3 

0 

3,7177 X 10~ 7 

6.1621 x lO' 5 

0 

0 

0 

M 

16.908 

16.908 

16.908 

8.836 

14.455 

16.649 

7.695 

9.228 

11.166 


16.908 

16.909 

16.974 

9.436 

14.450 

16.649 

7.695 

9.229 

11.529 

TT Btu 

-498.9 

-498.9 

-498.9 

12 183.9 

3820.2 

1900.3 

23 860.8 

14 872.7 

11 308.5 

H ’i^ 

-498.9 

-498.9 

-504.4 

11 292.9 

3826.2 

1901.3 

23 862.0 

14 876.6 

10 050.5 

H, 

-1.1596 x 10 8 

-1.1596 x 10 8 

-1.1596 x 10 8 

28.320 x 10 8 

8.8797 x 10 8 

4.4171 x 10 8 

55.462 X 10 8 

34.570 X 10® 

26.285 x 10 8 

kg 

-1.1596 x 10 8 

-1.1596 x 10 8 

-1.1724 x 10 8 

26.249 x 10 8 

8.8936 X 10 8 

4.4194 x 10 6 

55.465 x 10 6 

34.579 X 10 6 

23.361 x 10 6 

„ Btu 

3.9238 

3.1124 

2.5716 

6.9995 

4.0983 

3.2028 

8.1296 

5.6014 

4.3023 

lbm-°R 3 ' 9238 

3.1124 

2.5693 

6.6182 

4.0995 

3.2030 

8.1302 

5.6021 

4.1976 

O J 

16 417 

13 022 

10 760 

29 286 

17 147 

13 401 

34 014 

23 436 

18 001 

S ’ kg-°K 

16 417 

13 022 

10 750 

27 691 

17 152 

13 401 

34 017 

23 439 

17 563 
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contribution to existing knowledge. 
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contribution to existing knowledge. 
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